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phenomena, which specifically reflect the nature of the superconducting pairing state. 
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I. INTRODUCTION 

Understanding impurity effects was one of the early, clas- 
sical tasks for the theory of superconductivity f?.'?]. The An- 
derson theorem explained why non-magnetic impurity scat- 
terers do not destroy conventional, isotropic i-wave supercon- 
ductivity 1 3]. For unconventional superconductors, however, 
both non-magnetic and magnetic impurities are pairbreaking 
and the sensitivity to impurities and the concomitant suppres- 
sion of superconductivity may even serve as a fingerprint of 
an unconventional pairing state. It took three decades after 
the development of the BCS theory of superconductivity, be- 
fore experimental studies of impurity effects on microscopic 
length scales became possible - mostly due to the advance 
of scanning tunneling microscopy (STM). The need for ex- 
perimental probes with a high spatial resolution arose from 
the discovery of a continuously growing number of uncon- 
ventional singlet and even triplet superconductors; ruthenates 
and high-temperature cuprate superconductors belong to the 
most prominent examples. 

During the last decade it has furthermore become possible 
to manipulate surface structures on the atomic level. In par- 
ticular the achievements in the research group of D. Eigler led 
to the design of closed atomic arrangements of adatoms on 
metallic surfaces 0,0]. In these quantum corrals with circu- 
lar or elliptical shape electrons are confined to well defined 
geometries and STM techniques were applied to explore the 
structure of quantum mechanical wavefunctions and their in- 
terference patterns |5, 6]. Quantum mirage phenomena are 
observed, when an additional impurity atom is placed at one 
focus of an elliptic corral 0|. The terminology of quantum 
mirages is hereby defined by the observation that mirror im- 
ages of the local density of states pattern around the impurity 
at one focus point appear also at the second impurity free fo- 
cus. In subsequent studies even Kondo resonances could be 
detected at the impurity-free focus point, if a magnetic impu- 
rity adatom was placed in the other focus 18]. 

Similar mirage phenomena are expected for elliptic corrals 
on superconducting surfaces [9]. The local suppression of the 
superconducting order parameter around an impurity at a se- 
lected point in the corral should be transferred to image struc- 
tures due to standing wave patterns extending over the entire 
area of the elliptic corral. Superconducting quantum corrals 
therefore provide unique systems to study the combined ef- 



fects of impurity induced fingerprint structures of the pairing 
state and the mirage phenomena of interfering quantum me- 
chanical waves. With the assumption that a superconducting 
state is established in the bulk and on the surface of the sub- 
strate material we discuss in this paper the expected phenom- 
ena within the framework of BCS theory and the Bogoliubov- 
de Gennes equations. 

II. THE MODEL 

We consider the following mean-field Hamiltonian for an 
isotropic and an anisotropic singlet superconductor, respec- 
tively, with magnetic or non-magnetic impurities 

H = Y^U\Mfl (r)i/o¥a (r) + + H,„,p ( 1 ) 

with H{) — (— — /j) + ^(r). /J is the chemical potential, 
and y (r) is a hard-wall potential, which confines the electrons 
to the interior of an elliptic corral. The Schrodinger equation 
for particles in an elliptic geometry was previously solved an- 
alytically; the resulting eigenfunctions are a combination of 
Mathieu functions ccr, se, and modified Mathieu functions 
Ce,, lHollTlll: 

(P^_„^(e,1l) = Cer{QXn)Cer{T\Xn), (2) 

where 9 and r| are elliptical coordinates, r, n^^^^-^ enumerate 
the quantum numbers for the eigenstates. Their eigenenergies 

iHi 



are determined by the zeros of the modified Mathieu functions 

k'rfn ■ a denotes the length of the semimajor axis and e the 
eccentricity of the ellipse; m is the electron mass. 

For an isotropic s-wave superconductor with an attractive 
contact interaction //, is given by 

Hs = -j d^r [v|(r)vI(r)A(r)+/z.c.], (5) 

where A(r) = g{^\^ (r)) is the order parameter and g>Q 
the pairing interaction strength. Similarly, the mean-field 
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Hamiltonian for a singlet superconductor with a spatially ex- 
tended pairing interaction is described by 



d^r / dV [vI(r)vl(r')A(r,r') +/i.c.], (6) 



where the order parameter A (r,r') =g(r,r')(\|/^(r')\|/^(r)) de- 
pends on both coordinates r and r' of the electrons form- 
ing the Cooper pair Our ansatz for the pairing interaction is 
g(r,r') = g 5(|r — r'l —R) and assumes an attraction for elec- 
trons at a distance R\ this lengthscale R should be set by the 
size of a typical crystal lattice constant. Although g(r, r') is an 
isotropic real space interaction, it gives rise to an anisotropic 
local order parameter as we discuss below. The impurity po- 
tentials for magnetic (— ) and non-magnetic (+) impurities are 
described by 

= / d^r U{v) [v|(r)v^(r) ± v;(r)v^(r)] , (7) 



where t/(r) = L"L'o5(r — ri) with impurity positions r^. In 
this model for the magnetic scattering center the impurity 
spins are treated as classical spins, 5^1 which is a 

reasonable assumption e.g. for Mn and Gd adatom impurities 
on a niobium surface Ll4l1 . 

The mean-field Hamiltonian Eq. is diagonalized by the 
transformation 



¥a(r) = ^[M„(r)a„a + i^n(a)v„(r)a,t_J, 



(8) 



a„ 



d^r [M„(r)\(/o(r)-5g«(a)v„(r)i|/L<,(r)], (9) 



where a„o describes a superconducting quasiparticle with spin 
O and energy E„. The ansatz leads to the Bogoliubov-de 
Gennes equations 



[//o + f/(r)]M„(r)+ j dVA(r,r')v„(r')=£„M„(r), (10) 
dVA(r,r')M„(r')- [i/o±f/(r)]v„(r) -£„v„(r), (11) 



where the +1- sign in Eq. ( II 1> corresponds to non- 
magnetic/magnetic impurities, respectively. These equa- 
tions are simplified in the case of an isotropic superconduc- 
tor, where the pairing-interaction is pointlike, and therefore 
A(r,r') ~ 5(r — r') A(r). The order parameter A(r) is deter- 
mined self-consistently from the condition 



A(r)=g^«„(r)v„(r) [l-2/(£„)], 



(12) 



where /(£) denotes the Fermi function. For an anisotropic 
superconductor with a finite range pairing interaction we have 
instead 

A(r,r') =g(r,r')I^{M„(r)v„(r')[l -/(£„)] (13) 

n 

(r')v„(r)/(£„)}- 

In order to solve the Bogoliubov-de Gennes equations ([lO} 
and (lll> we expand u„ (r) and v„ (r) in terms of free electron 



eigenfunctions inside the elliptic corral 111 
«"(«■) = (P/t(r), 

k 



(14) 
(15) 



where k enumerates all the eigenstates given in Eqs. (|2j and 
(|3}- For the ground-state energies of the isotropic (Eq. il6i ) 
and the anisotropic (Eq. ilH ) superconductor we find 

^s- /d^r [y„(r)//oy„(r)-A(r)M„(r)y„(r)], (16) 



2v„(r)//ov,,(r) 



dVA(r,r')(M„(r)v„(r')+M„(r')v„(r)) 



(17) 



The local density of states (LDOS) is obtained from 

A^(r, CO) = ^ [«'(r)5((0 - E„) + v2(r)5(co + £„)] , (18) 

n 

and the total density of states (DOS) is N{(x)) ^ J d^rN{r,&). 

III. RESULTS 

For the elliptic corral we choose its eccentricity e = 0.5 and 
the length of the semi-major axis a = 150 A; these two param- 
eters determine the single-particle eigenenergies and -states. 
With an approximate electronic density of one electron per 
6Ax6A square the total number of electrons in the corral is 

— 1726. With this choice the 863"' eigenstate at the Fermi 
energy has a large probability density at the foci, which 
is a favorable situation for the occurrence of mirage phenom- 
ena. For an isotropic superconductor the mirage phenomenon 
is observed, whenever there is one of these non-zero proba- 
bility density eigenstates lying inside the energy gap. For our 
choice of the pairing interaction strength g = 1.0 x lO^^Ef 
the energy gap Aq stretches over about 50 eigenstate ener- 
gies, where several of these states have the required property. 
Hence, the existence of mirage phenomena for an isotropic 
superconducting surface does not depend sensitively on the 
position of the Fermi energy, in contrast to a metallic, normal 
conducting surface fioll and also in contrast to an anisotropic 
superconductor, as we discuss below. 

A. Isotropic .s-wave superconductor 

We start our analysis with a corral on the surface of an 
isotropic i-wave superconductor In the presence of non- 
magnetic impurities the system remains time-reversal invari- 
ant. These impurities do not break Cooper pairs, which are 
built from Kramers-degenerate electronic states [3]. Conse- 
quently, no additional states are created in the renormalized 
energy gap. Magnetic impurities on the contrary, are pair- 
breaking defects, because the spin-up electron of the Cooper 
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FIG. 1: Local density of states (LDOS) for an isotropic supercon- 
ductor with a magnetic impurity at the right focus. Dashed and 
dash-dotted lines show the LDOS at the impurity site (0.5a, 0), dot- 
ted and dash-double-dotted lines show the LDOS at the mirror point 
(—0.5a, 0) . The arrows ascribe the spin projection to the particle and 
the hole component of the bound state. 



FIG. 2: Bound-state energy for an elliptic quantum corral (triangles) 
with a magnetic impurity at the right focus (0.5a, 0) compared with 
the bound-state energy for a magnetic impurity on an open surface 
of a i-wave superconductor (solid curve); the arrows mark the spin- 
direction of the particle component of the quasiparticle bound state. 
Filled/unfilled triangles correspond to a spin-down/spin-up quasipar- 
ticle, respectively; Uq/ef = 0.03(/o/Ao- 



pair is repelled (Uq > 0) and the spin-down electron is at- 
tracted by the impurity, leading to a quasiparticle bound state 
with a spin-down particle peak and a spin-up hole peak in the 
LDOS shown in Fig. ^(dashed and dotted curves), where the 
impurity is located at the right focus of the ellipse. The LDOS 
at the impurity site vanishes at the gap edge energies, while the 
LDOS at the mirror point (dotted line) is finite at these ener- 
gies. The total spectral weight of both peaks, i.e. the spatially 
and energy (between the gap edge energies — Aq and Aq) inte- 
grated LDOS is unity, characterizing one localized quasiparti- 
cle. Hence, the bound state corresponding to the dashed line is 
a local spin-down quasiparticle state (see also Fig. |3 unfilled 
triangles), where the spin projection of the particle component 
is antiparallel and that of the hole component is parallel to the 
impurity spin. This is the same phenomenon known for the ex- 
tended s-wave superconductor without confining boundaries 
fisll . The dash-dotted and dash-double-dotted curves in Fig. 
n] show the LDOS for a stronger impurity potential. The dif- 
ference now is that the spins of the particle and the hole peak 
have changed, i.e. we find a spin-up quasiparticle (compare 
with Fig. |2j filled triangles). This change is due to the fact 
that both peaks cross the Fermi energy at a critical impurity 
potential strength Uc- 

In the corral geometry we observe the two-peak structure in 
the LDOS not only at the impurity site e.g. at one focus point, 
but also - with some attenuation - at the impurity free focus 
(dotted and dash-double-dotted lines) |9]. The mirage effect is 
only strong at one of the two peaks at the mirror point LDOS. 
There the particle peak carries very little spectral weight for 
Uo < Uc (see Fig. ^ dotted line), and so does the hole peak for 
Uo > Uc- But with increasing Uq the particle and hole peaks 
at the mirror point are both growing. We emphasize that the 



LDOS peaks at the mirror point are at the same resonance en- 
ergies as in the LDOS at the impurity site, because they belong 
to the same bound state with energy coq. This is the quantum 
mirage effect in the LDOS for an elliptic corral built on the 
surface of an isotropic i-wave superconductor As expected, 
we find this effect to be very robust against a change in the 
Fermi energy and thus the electron density. Moreover this mi- 
rage effect remains almost unaffected by an additional pertur- 
bation in form of a second magnetic impurity placed anywhere 
else inside the corral. 

Since the bound-state peaks at energies co = icoo, split 
off from the gap edges, they move with increasing potential 
strength Uo symmetrically towards the Fermi energy (see Fig. 
|2}. coo therefore decreases until it reaches the chemical poten- 
tial for a critical potential Uc, which depends on the impurity 
position ro. At Uo = Uc the bound state becomes a zero energy 
state, and the ground state of the superconductor becomes un- 
stable. The existence of this critical point, which signals a 
ground-state level crossing transition, was first pointed out by 
Sakurai for i-wave superconductors with magnetic impurities 
lITsll . At Uc the system undergoes a first order phase tran- 
sition from a spin zero to a spin — j ground state fl?, Tsf] . 
While the condensate ground state for Uo < Uc consists only 
of Cooper pairs, the new ground state for Uo > Uc contains an 
additional spin-down electron (or up-spin hole). The impurity 
induced bound state is now a spin-up quasiparticle (Fig. Q] 
dash-dotted and dash-double-dotted curves and Fig. |2j filled 
triangles). A further increase of Uo in this regime (Uo > Uc) is 
accompanied by a spatial shift of the LDOS at the resonance 
energies away from the impurity position. This spatial redis- 
tribution of the spectral weight leads to decreasing heights of 
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the bound-state peaks in the LDOS Fig. [2 at the impurity po- 
sition (right focus) and an accompanying increase of the peak 
heights at the impurity free focus. Note that the total spectral 
weight of the bound-state peaks always remains unity. Hence 
we observe a transition from a mirage-effect for small Uo to 
an anti-mirage-effect for large Uq. 

We compare our self-consistent solutions of the 
Bogoliubov-de Gennes equation for a corral with the re- 
sults for an open surface using the following formula for the 
bound-state energies 



coo = Ao 



\l-i%NFfU^\ 
1 + (nNpfU^ 



(19) 



where Nf is the normal state DOS at the Fermi energy. Eq. 
( I19> is obtained from the non-self-consistent T-matrix formal- 
ism, which proved to be sufficient and even quantitatively ac- 
curate for the description of magnetic impurity induced bound 
states in i-wave superconductors without confining walls f3. 
In Fig. (|2ji we use Eq. ( I19l l with kNf Nf (solid curve), 
where Nf desribes the mean density of states averaged over 
an energy window ± Aq in the normal state. Fig. |2]shows 
that the bound-state energies in the elliptic corral are well de- 
scribed by Eq. (19). Comparing the self-consistent and non- 
self-consistent results for our corral, we find indeed only small 
differences: neither the quasiparticle energies are shifted no- 
tably nor the heights of the quasiparticle peaks are modified 
significantly. 

Another interesting quantity is the renormalized local or- 
der parameter A(r) given by Eq. il2\ . We observe a sup- 
pression of A(r) at the impurity site, where the lengthscale 
of the suppression and the wavelength of the order parameter 
oscillations is given approximately by the Fermi wavelength 
Xf = Inh/ ^ImEp ~ 15 A, but seems to be independent of 
the size of the energy gap Aq. Remarkably, the local order pa- 
rameter at the impurity position changes discontinuously from 



a positive to a negative value, when Uq crosses the critical 
potential strength Uc lll3l [ItIi . For a non-magnetic impurity 
the reduction of A(r) in the vicinity of the impurity is due to 
the repulsive potential for both spin directions resulting in a 
reduced electronic density at and near the impurity site and 
therefore to a reduced local pairing amplitude. In Fig. |3]a 
magnetic impurity is placed at the right focus point. All one- 
particle eigenstates with a large probability density at the right 
focus have a large probability density at the left focus as well. 
The difference plot of the order parameter Fig. |3]strongly re- 
sembles the 863''' eigenstate, implying that A(r) is suppressed 
not only at the impurity site (right focus), but almost symmet- 
rically at the impurity free (left) focus, which we therefore 
call an order parameter mirage effect. We note that the order 
parameter is negative at the impurity site in Fig. |3] 

If there are magnetic impurity spins in both foci, their an- 
tiparallel alignment leads to a lower ground-state energy Eq. 
( I16> . Moreover we observe that the order parameter suppres- 
sion is about six times stronger for a parallel than for an an- 
tiparallel impurity spin alignment. A parallel alignment leads 
to a four-peak structure in the LDOS, because there are two 
similar states localized at each impurity, which are hybridiz- 
in g an d therefore splitting into bonding and antibonding states 
B flql . If the impurity spins are aligned antiparallel, we find 
only two peaks inside the energy gap. The intensity of the 
hybridization depends in general on the distance between the 
two impurities; for a corral geometry this dependence is quite 
complex. 



B. Anisotropic superconductor 

Now we assume the pairing interaction ^(r,r') to be at- 
tractive in a distance of a typical crystal lattice constant R ^ 
0.04a = 6A. In this case the order parameter A(r,r') depends 
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FIG. 4: Local order parameter A(r,r + r') in units of Zf for the anisotropic superconductor for a fixed r plotted as a function of r' = {x' ,y') 
with y,y e [—2R,2R]; the dashed lines mark the order parameter nodes; r = (0,0) in the left panel and r = (0.44a, 0) in the right panel. 



on the position inside the corral as well as on the relative coor- 
dinate. Although the interaction is still isotropic, the order pa- 
rameter is now anisotropic. In Fig. |3(left panel) we show the 
order parameter in the center of the ellipse, where we identify 
a d-wave like structure with sign changes and nodes at angles 
near 45° with respect to the semimajor axis. At the center 
the order parameter is least affected by the confining corral, 
but except for this point of highest symmetry, the structure 
of the local order parameter is distorted in a complex way. 
An example for a remarkably complicated symmetry and lo- 
cal structure of A(r,r') is shown in the right panel of Fig. |3 
The spatially extended pairing interaction leads to a V-shape- 
resembling DOS (see Fig. |5] solid curve), with no sharp gap 
edges. Furthermore, the DOS has no particle-hole symmetry, 
because of the asymmetric distribution of eigenenergies for 
the corral eigenfunctions. 

While non-magnetic impurities do not give rise to local- 
ized states in a ^-wave superconductor, they do have a pair- 
breaking effect in an anisotropic superconductor |18]. For 
a non-magnetic impurity placed at different positions inside 
the corral we calculate the local order parameter, the DOS, 
and the LDOS. In the DOS we always find bound-state peaks, 
which move towards the Fermi energy as Uq increases, but 
never reach . Hence, we do not observe a zero-energy-peak 
in agreement with the results for an open surface in a particle- 
hole asymmetric situation 1 13i 191 . The DOS is modified over 
the entire energy range as shown in Fig. |5l For weak impurity 
potentials we observe a two-peak structure in the DOS near 
the Fermi energy, but after exceeding some critical potential 
Uc, which depends again on the impurity position ro, each of 
the two peaks starts to split into two, so that a four-peak struc- 
ture emerges near co = Ef (see Fig. |5jl. At the resonance en- 
ergies for small Uq (two-peak regime) we find, that the bound 
state has a large LDOS in the vicinity of the impurity. By 
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FIG. 5: Density of states (DOS) for the anisotropic superconductor 
with a non-magnetic impurity at (0.44a, 0). 



increasing Uq the original bound state acquires admixtures of 
the 862"'' and ^66"" state f20^. Above the critical potential Uc, 
we identify the 866''' state either at both outer or both inner 
peaks of the four-peak low-energy DOS, while the LDOS at 
the two other peaks consists of the small Uq bound state mixed 
with a contribution of the 862"'' state. Hence the splitting of 
the bound state energy corresponds to a "demixing" of the 
contributing states, which form the bound state. We observed 
this behavior at different impurity positions on the semimajor 
axis, but the occurrence of a four-peak structure nevertheless 
depends sensitively on the Fermi energy. Slight variations of 
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the electronic densities can lead to two-peak structures in the 
DOS, irrespective of the impurity potential strength. 

In order to show a few of a large variety of interesting and 
astonishing effects we calculate the LDOS in the presence 
of a non-magnetic impurity in particular at the right focus 
(0.5fl,0), at (0.44a, 0), and at the center (0,0). Interestingly, 
the important states (862"'' and 866''') mentioned above are 
the only contributions to the bound state for an impurity in 
a focus point. These states have a zero LDOS on the entire 
semimajor axis, but several maxima about a lattice constant 
away from the focus. This bound state shows an almost per- 
fect reflection symmetry with respect to the semiminor axis. 
Shifting the impurity a distance 1.5R away from the right fo- 
cus to (0.44fl, 0), we observe that the LDOS at the lower reso- 
nance energies (a^ = ±2.8 x lO^^Ef has a strong contribution 
of the 863"' state and, as discussed before, a weaker contribu- 
tion of the 862"'' state. Fig.|6]shows the LDOS at the positive 
resonance energy co^ . While the LDOS at the impurity site it- 
self is zero, we observe a maximum in the LDOS at the mirror 
point. The impurity induces asymmetrically spectral weight 
around both foci, whereas the height is larger at the left than 
at the right focus, which is close to the impurity. Thus, we ob- 
serve a surprising, strong anti-mirage effect in the LDOS at all 
resonance energies with a non-magnetic impurity at (0.44a, 0) 
independent of the potential strength Uo- 

In Fig. 0we show the LDOS with an impurity located at 
the center of the ellipse at the energy of the higher positive 
resonance energy co^ — 7.6 x lO^^e^ . We identify the 862""^ 
eigenstate again, which determines the region near the bound- 
ary of the ellipse and consists of rings of peaks aligned along 
four ellipses. The structure of the LDOS in the vicinity of the 
impurity is in fact more influenced by the symmetry of the lo- 
cal order parameter than the geometry of the problem. E. g. 
in a range of a few lattice constants around the center there 
is no LDOS along the diagonal directions, and locally there 
is almost a fourfold symmetry as for the absolute value of the 



order parameter in Fig. |4](left panel) i2lll . 

In the vicinity of the impurity site the order parameter is 
reduced strongly by a non-magnetic impurity. This becomes 
obvious by comparing the right panels of Figs. |4]and|8] where 
the impurity is located at (0.44fl,0). Moreover we observe 
here and quite generally a modificaton of the spatial structure 
and the symmetry of the local order parameter, which depends 
sensitively on the potential strength Uq. Interestingly the new, 
impurity induced spatial structure resembles more a d-wave- 
like symmetry than the original structure. Fig. |8] clearly in- 
dicates, that this t/-wave-like order parameter structure is pro- 
jected to the impurity free mirror point. Instead of a reduction, 
the order parameter is actually enhanced at the mirror point. 
As in the LDOS we therefore also observe an anti-mirage ef- 
fect for the order parameter. The order parameter at (0.44fl, 0) 
is changed to a structure similar as in Fig. |8](right panel), if an 
impurity is located at the center (!) of the ellipse. In general 
we find, that the magnitude of the order parameter is not only 
reduced over the entire area of the corral, but also modified 
with respect to its local spatial symmetry. 

Regarding the connection between the spatial structure of 
the LDOS in the vicinity of the impurity and the order param- 
eter, we encounter three different situations, (i) the LDOS re- 
sembles the order parameter in the absence of an impurity, (ii) 
the LDOS resembles the order parameter in the presence of 
an impurity, and (iii) there is no connection at all between the 
two. The latter is observed for example for a focus point impu- 
rity. We find the behavior (i) only for an impurity at the center 
of the corral, for large impurity potentials, and at positive res- 
onance energies. For smaller impurity potentials we are back 
to case (ii), where again the LDOS shows the structure of the 
order parameter with an impurity only at positive resonance 
energies. This second case seems to reflect the more generic 
situation. Hence, the LDOS in the presence of a non-magnetic 
impurity inside the corral can directly reveal the structure and 
the local symmetry of the order parameter. 
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FIG. 8: Local order parameter A(r,r + r') for a fixed r with a non-magnetic impurity at (0.44fl,0); Uq/Ef — 2.8 x 10^^; r — (— 0.44a, 0) in 
tlie left panel and r = (0.44a, 0) in tlie riglit panel. See also the caption of Fig. 4. 



IV. CONCLUSION 

We have observed a rich variety of mirage- and anti-mirage- 
phenomena not only in the LDOS, but also for the local order 
parameter structure. The induced patterns in both quantities in 
the vicinity of an impurity lead to mirror images in impurity- 
free regions, which are characteristic for both, the structure of 
the electronic wavefunctions in the elliptic geometry and the 
nature of the superconducting state. With the continuing ad- 
vances of STM techniques and the already achieved capabil- 
ities for the design of atomic corral arrangements on metallic 
surfaces the realization of quantum corrals on superconduct- 



ing surfaces indeed appears to be possible. From our model 
analysis we are led to expect intriguingly rich structures in 
the LDOS patterns, which by themselves contain information 
about the nature of the superconducting state. 
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